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Abstract. We consider the dephasing of two internal states |0) and |1) of a trapped 
impurity atom, a so-called atomic quantum dot (AQD), where only state |1) couples 
to a Bose-Einstein condensate (BEC). A direct relation between the dephasing of 
the internal states of the AQD and the temporal phase fluctuations of the BEC 
is established. Based on this relation we suggest a scheme to probe BEC phase 
fluctuations nondestructively via dephasing measurements of the AQD. In particular, 
the scheme allows to trace the dependence of the phase fluctuations on the trapping 
geometry of the BEC. 
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1. Introduction 

The coherence properties of Bose-Einstein condensates (BECs) have attracted 
considerable theoretical and experimental interest since the first experimental realisation 
of BECs in trapped ultracold clouds of alkali atoms [H |2] . Part of this interest is due to 
the importance of coherence effects for the conceptual understanding of BECs and their 
use as a source of coherent matter waves. In particular, the absence of spatial coherence 
in low dimensional BECs, which exhibit strong spatial and temporal phase fluctuations, 
has been investigated theoretically |S1 IS Ej and demonstrated for one dimensional 
condensates [HI Ej. Moreover, temporal first and second order phase coherence of an 
atom laser beam extracted from a BEC has been observed jHl IHI- However, to our 
knowledge, temporal phase fluctuations in a BEC have not yet been measured directly, 
despite the fact that they ultimately limit the coherence time of an atom laser beam. 

In this paper, we propose a scheme to measure temporal phase fluctuations of 
the BEC based on a single trapped impurity atom coupled to the BEC, hereafter 
called an atomic quantum dot (AQD) jTHI. More specifically, we consider an AQD 
with two internal states |0) and where we assume for simplicity that only state |1) 
undergoes collisional (s-wave scattering) interactions with the BEC. This setup could 
be implemented using spin- dependent optical potentials [TT], where the impurity atom, 
trapped separately ^21 ^1 BEC atoms correspond to different internal 

atomic states. The collisional properties of the AQD can, to a good approximation, 
be engineered by means of optical Feshbach resonances ^21, which allow to change 
the atomic scattering length locally over a wide range and are available even when no 
magnetic Feshbach resonance exists. 

By identifying the combined system of the AQD and the BEC with an exactly 
solvable independent boson model [Tn|II7l|TH], we show that the dephasing of the internal 
states due to the asymmetric interaction with the BEC is directly related to the temporal 
phase fluctuations. This dephasing can be detected under reasonable experimental 
conditions, for example, in a Ramsey type experiment and hence it is possible 
to use the AQD to probe BEC phase fluctuations. Since the phase fluctuations depend 
strongly on the temperature and the density of states of the BEC, determined by the 
trapping geometry, the proposed scheme allows us to measure the BEC temperature and 
to observe the crossover between different effective BEC dimensions, notably transitions 
from 3D to lower dimensions. Our scheme is nondestructive and hence it is possible, 
in principle, to investigate the dependence of phase fluctuations on the BEC dimension 
and temperature for a single copy of a BEC. 

Probing a BEC with an AQD was proposed recently in ^U] and in a different 
context in [201 ■ However, in jT^ two states corresponding to the presence of a single 
atom in the trap or its absence were considered, as opposed to internal atom states. 
More importantly, in addition to the collisional interactions the impurity atom was 
coupled to the BEC via a Raman transition, allowing the realisation of an independent 
boson model with tunable coupling. In particular, Recati et al [TU] proposed to measure 
the Luttinger liquid parameter K by observing the dynamics of the independent boson 
model for different coupling strengths. 
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Other schemes based on interactions between impurity atoms and BECs have been 
proposed recently. In [21], for example, an impurity was used to implement a single 
atom transistor, whereas in [22] a quantum gate exploiting phonon-mediated interactions 
between two impurities was investigated. Single atom cooling in a BEC was considered 
in j2S] where some aspects of dephasing of a two state system (qubit) due to BEC 
fluctuations were addressed. Whereas their treatment focused on a three-dimensional 
BEC and was based on a master equation approach, in our work we use an analytical 
approach similar to dephasing calculations for semiconductor quantum dots [211 123 • 

The paper is organized as follows. In section 2 we introduce our model for the 
AQD coupled to the BEC and identify it with an exactly solvable independent boson 
model, which enables us to establish a relation between the temporal phase fluctuations 
of the BEC and the dephasing of the AQD. We first consider an effective D-dimensional 
BEC, which is strongly confined in (3 — D) directions and, as an ideal case, assumed 
to be homogeneous along the loosely confined directions. Later we discuss corrections 
to this zero potential approximation due to a shallow trap potential. In section 3 we 
show by means of a concrete example how dephasing depends on the interaction time 
r between the AQD and the BEC, the condensate temperature T and the effective 
condensate dimension D. Section 4 addresses corrections to our results due to an 
imperfect measurement of the AQD state. We conclude in section 5. 

2. The model 

We consider the system of an AQD coupled to a BEC in thermal equilibrium with 
temperature T. The BEC is confined in a harmonic trapping potential 



with r = (x, y, z) the position vector and m the mass of the condensate atoms. The 
trap frequencies oj^, ojy and Uz assume the value(s) uj±_ or where ci;j_ 3> a; is the trap 
frequency in the strongly confined (transverse) direction(s) and uj the trap frequency 
in the loosely confined direction(s). Depending on the number of strongly confined 
directions - none, one or two - the BEC is either 3D or assumes effective 2D or ID 
character, provided that thermal excitations are suppressed, i.e. ksT -C huj±, and that 
the interaction energy of the weakly interacting BEC does not exceed the transverse 
energy, i.e. mc^ -C huj^, with fc^ Boltzmann's constant and c the speed of sound. In 
directions that are loosly confined the extension of the trap potential is assumed to 
be much larger than the length scale a set by the AQD size, so that it is justified to 
approximate the potential by zero, as discussed at the end of this section. We consider 
the case where the spectrum of the BEC excitations, as compared to the impurity 
spectrum, is practically continuous. 

The AQD consists of an impurity atom in the ground state of a harmonic trap 
potential centered at ro and is described by the wave function ifjair—Yo). We assume that 
ipa[T — VQ) takes the form of the BEC density profile in the strongly confined direction(s) 
and has ground state size a in the loosely confined direction(s). For example in case of 




(1) 
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a ID condensate 

V'.(r-ro)oc^=exp =— j - (^^J - (^=— j , (2) 

with the harmonic oscillator length a± <^ a. We further assume that the impurity atom 
has two internal states |0) and |1), and undergoes s-wave scattering interactions with 
the BEC atoms only in state |1). 



2.1. Dephasing in the zero potential approximation 
The total Hamiltonian of the system can be written as 

Htot = Ha + Hb + Hi , (3) 

where Ha = with level splitting hVL., is the Hamiltonian for the AQD, Hb is the 

Hamiltonian for the BEC, and Hj describes the interaction between the AQD and the 
BEC. Provided that a// ^ 1, with / the average interparticle distance in the BEC, we 
can represent the D-dimensional (quasi) condensate in terms of the phase operator (/>(x) 
and number density operator ntot(x) = tiq + n(x). Here, no = is the equilibrium 
density of the BEC, n(x) is the density fluctuation operator and x a D-dimensional 
vector describing the position in the loosely confined direction(s). We describe the 
dynamics of the BEC by a low-energy effective Hamiltonian [2^1 HZ] 

HB = \j d^x(^no(V0)2(x) + gh\^)j , (4) 

where the interaction coupling constant is defined hy g = mc^/riQ |. We shall see 
that the use of this model is fully justified, even though it exhibits only Bogoliubov 
excitations with a linear dispersion relation ujk = ck, henceforth called phonons. The 
canonically conjugate field operators 0(x) and n(x) can be expanded as plane waves 

0(x) = Y: A,\b]^ exp [ik . x] + h.c. ) , (5) 



n X 



exp [ik ■ x] - h.c. ) , (6) 



where and b]^ are bosonic phonon creation and annihilation operators and is the 



sample size. With the amplitudes = yhuk/g the Hamiltonian takes the familiar 
form 

HB = j:^Mblb], + h. (7) 
k ^ 

The coupling between AQD and the BEC occurs in the form of a density-density 
interaction 

Hi = k\1){1\ /d^x|^.(x-Xo)pntoi(x), (8) 



where ipa{'^ — xq) is the AQD wave function integrated over the transverse direction(s), 
and K the coupling constant, which can be positive or negative. To avoid notable 

I We neglect the weak dependence of c on the dimensionahty of the BEC |2H| for simpHcity. 
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deviations of the BEC density from hq in the vicinity of the AQD we require |k| ~ 
[2niEni- The interaction Hamiltonian Hj acts on the relative phase of the two internal 
states, but does not change their population. This effect is customarily called pure 
dephasing. 

The total Hamiltonian can be identified with an independent boson model, which 
is known to have an exact analytic solution ^HIEIEI- Inserting the explicit expression 
(jni for n(x.) into equation (jH} we can rewrite the total Hamiltonian as 

Htot = i^no + hn)\i){M + E {9 A + 4^k) |iXi| + E ^^^(^L^k + h . (9) 

k k ^ 

where khq is the mean field shift. The coupling coefficients g-^^, which contain the specific 
characteristics of the system, are given by 

^k = -;7§n^^4' (10) 
with the Fourier transform of the AQD density 

/k = y'd^x|^,(x-Xo)|2exp[-ik.x]. (11) 

As can be seen from expression © the effect of the AQD in state |1) is to give rise to 
a displaced harmonic oscillator Hamiltonian for each phonon mode. This problem is 
identical with the problem of a charged harmonic oscillator in a uniform electric field 
and can be solved by shifting each phonon mode to its new equilibrium position ^H] ■ As 
a consequence the AQD does not change the phonon frequencies and hence the phase 
fluctuations are the same in the presence of the AQD as in its absence. 

The state of the system is described by the density matrix ptot(t). We assume that 
the AQD is in state |0) for t < and in a superposition of |0) and |1) at t = 0, which 
can be achieved, for example, by applying a short laser pulse. Hence the density matrix 
Ptot{t) at time t = is of the form 

PtotiO) = p(0) (E)PB, PB = ^ exp[-HB/{kBT)] , (12) 

^B 

where p(0) is the density matrix for the AQD and pb the density matrix of the BEC 
in thermal equilibrium, with Zb the BEC partition function. After a change to the 
interaction picture the time evolution operator of the total system takes the form ^7] 



U{t) = exp 



k 



(13) 



with = g]^{l — exp[iukt])/{huJk)- The reduced density matrix of the AQD can be 
determined from the relation p(t) = Tt b{U (t) ptot{0)U~^ (t)} , where Tr^ is the trace over 
the BEC. The coherence properties of the AQD are governed by the off-diagonal matrix 
elements 

Pio(t)=pSiW=Pio(0)e-^W, (14) 
with the dephasing function 



7(t) = — In <^ exp 



E K^L - Pik 

k 



(15) 
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where angular brackets denote the expectation value with respect to the thermal 
distribution ps- 

The dephasing function 7(t) can be expressed in terms of the phase operator in the 
interaction picture 0(x, t) = exp[iifBt/^]0(x) exp[— iif^t/Zi]. We introduce the coarse- 
grained phase operator averaged over the AQD size a 



^(xo,t) = / d^x|V'.(x-Xo)|'0(x,t) 



and the phase difference 

as 

P-7W - 



(16) 



6o-(xo, t) = 0(j(xo, t) — 0o-(xo, 0) to rewrite the phase coherence 



exp 



exp 



i-50<^(xo,t) 
9 

1 
2 



(17) 



K\2 



9 



(xo,t) 



where the second equality can be proven by direct expansion of the exponentials jl8j . 
Thus we have established a direct relation between the dephasing of the AQD and the 
temporal phase fluctuations of the BEC averaged over the AQD size a. Moreover, the 
phase coherence ffTHj) is closely related to the correlation function (\E'^(x, t)\E'(x', t')), with 
\E'(x, t) the bosonic field operator describing the BEC in the Heisenberg representation. 
In the long-wave approximation, where the density fiuctuations are highly suppressed, 
we have [23 123 

1 



(^t(xo,t)^(xo,0)) ^ no exp 



(50)1x0, t) 



(19) 



Therefore expression (fTHj) can be interpreted as a smoothed out temporal correlation 
function of the BEC field operator \E'(x, t) provided that n = g. 

To further investigate the effect of the BEC phase fluctuations on the AQD 
we require an explicit expression for the dephasing function 7(t) depending on the 
parameters of our system. For this purpose we express the dephasing function 7(t) in 
terms of the coupling coefficients g^ [T7j 

huk \ 1 — cos(ti;A,i) 



7(i) = El^?kl'coth(- 
k 



2kBT^ {hujky 



(20) 



and take the thermodynamic limit of 7(t), which amounts to the replacement 

/o°° dk g{k), with the density of states g{k) = Sd L^k^-\ Sd = D/[2^7r^/^T{D/2 + 1)] 

and r(x) the gamma function. Substituting expression (fTUI) for g^^ we flnd 



(21) 



-2kBT^ huok 

The integral is well deflned due to the factor = exp[— o"^A;^/2], which provides 

a natural upper cut-off. Since for typical experimental parameters ^/l ~ 1, with 
~ h/{mc) the healing length, the condition a /I ^ 1 implies that a ^ ^. Thus 
the upper cut-off at A; ~ l/cr is still in the phonon regime, which justifles the use of the 
effective Hamiltonian (j^. 
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2.2. Corrections due to a shallow potential 

We now give the conditions under which the effects of a shallow trap on result ()2H) are 
negligible in the thermodynamic limit. This limit is taken such that uj oc 1/-R, with R 
the radius of the BEC, to ensure that the BEC density at the center of the trap remains 
finite. To determine the trap-induced corrections to the density of states and the phonon 
wave functions, which affect /j^ defined by (fTT|) . we use a semi-classical approach based 
on the classical Hamiltonian [3T^ 

/f(p,x) = c(x)|p| + V(x), (22) 

with c(x) = c(l — x^/i?^)^/^ the position dependent speed of sound, p the phonon 
momentum and V(x) = muo'^y? /2 the shallow trap potential. The range of phonon 
energies e relevant for the AQD dephasing is huj/2 < e < hc/a, where the lower bound 
goes to zero in the thermodynamic limit. 

Given the semi-classical phonon wave functions, we find that corrections to /j^ are 

negligible if a <^ Lg, with = \J'2,e/ {mu'^) the classical harmonic oscillator amplitude. 
For the density of states g{e) we use the semi-classical expression 

9{e) = I d^x J d^p 6{e - H{p, x)) (23) 

to obtain g{e) oc e^^^ [1 + 0{Ls/R)]. Thus in the regime where a <^ <^ R and 
after the substitutions L ^ and huk — > e we have 

oc-/ d££^-V.I'coth(— -) (24) 

g Jo ^ikbi ' s 

which up to numerical constants is identical to expression (|2H). The conditions on L^, 
together with the bounds of the phonon spectrum, imply that the dephasing in a shallow 
trap does not differ from the homogeneous case if ^ -C a -C a^j, with a^j = \Jh/{muj) 
the harmonic oscillator length. 

3. Application 

In this section we discuss by means of a concrete example § how the AQD dephasing 
depends on the interaction time r, the condensate temperature T and the effective 
condensate dimension -D ||. To find the phase coherence e"'^'-^'' we have evaluated the 
integral ()21|) numerically (shown in figures 1, 2 and 3) and analytically in the high 
temperature regime fc^T ^ hc/a. 

Figure ^ shows the phase coherence e"'*''^'^-* as function of the interaction time r. 
The evolution of the phase coherence is split into two regimes separated by the typical 
dephasing time cj/c ~ 10~^s, which can be viewed as the time it takes a phonon to 

§ For the numerical values of the system parameters see the caption of figure ^ 

II A discussion of a similar model in the context of quantum information processing can be found in 
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Figure 1. The phase coherence e"^^'^-' as a function of the interaction time t. In 
the regime t ^ a /c the phase coherence approaches a non-zero value in 3D (sohd), 
faUs off as a power law in 2D (dashed) and shows exponential decay in ID (dotted). 
The system parameters were set to m = 10~^^kg, Z = 5 x 10~^m, c = 10~'^ms~"'^, 
T = 2 X lO-'^K, a = lO-^m and n^g. 



pass through the AQD. It follows from the analytic calculations that the asymptotic 
behaviour of e"'^'-^'' for times r ^ a/c is 

'K.\2 mckBT 

'9' 



-7(r) 



Ct exp 



for D 



exp 



K\2 mkBT 



■(-) 



a 



for D 



for D 



(25) 



where u = {k/ gYrnksT / and Ct, are temperature dependent constants. 
Thus the phase coherence tends asymptotically to a non-zero value in 3D, falls off as a 
power law in 2D and shows exponential decay in ID. This result reflects the fact that 
the physics of ID and 2D condensates differs significantly from that of a 3D condensate. 
Given the relation between the phase coherence and the BEC correlations, discussed in 
section 2.1, we note that the results for the phase coherence are consistent with findings 
for temporal and spatial coherence properties of the BEC [221 HZ] , where in the latter 
case the distance |x| has to be identified with cr. 

In particular in 3D the dephasing of the AQD is incomplete under reasonable 
experimental conditions due to the suppressed influence of low frequency fluctuations, 
which has been noted in the context of quantum information processing f7| 1^ . 
However, the residual coherence of the AQD will eventually disappear due to processes 
as, for example, inelastic phonon scattering [21], which are not taken into account in 
our model. We also point out that, in contrast to the remnant coherence of the AQD, 
the typical dephasing time does not depend on the coupling constant k. 

The fact that phase fluctuations depend strongly on the density of states, or 
equivalently the dimension of the BEC, allows us to observe the crossover between 
different effective dimensions, especially transitions from 3D to lower dimensions. 
Figure |21 shows the phase coherence e"'*''^'^'' as a function of the dimension D for the 
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Figure 2. The phase coherence e"''^'^-' as a function of the effective condensate 
dimension D for the interaction times ti = cr/c (dotted), T2 — 2a/c (dashed) and 
T3 = 10 cr/c (solid). The phase coherence drops significantly while the BEC excitations 
in the strongly confined directions are frozen out. The plot was produced with the 
same set of parameters as in figure ^ 



interaction times ti = a / c, T2 = 2 cr/c and = 10 cr/c. The phase coherence e"''''-'^^ drops 
significantly as the BEC excitations in the strongly confined direction(s) are frozen out. 
Therefore we expect that the change of phase fiuctuations depending on the effective 
dimension should be experimentally observable. 

In addition, the AQD can be used to measure the BEC temperature since the 
dephasing function 7(r) is approximately proportional to the temperature T according 
to equation ()25|1 . Figure El shows the relation between the phase coherence e"^^"^-* and 
the BEC temperature T in all three dimensions for the interaction times ri = cr/c, 
T2 = 2a/c and = 10 cr/c. In ID and 2D the interaction time r can be chosen to 
assure that the AQD dephasing changes significantly with temperature, whereas in 3D 
the interaction time r has little infiuence on the thermal sensitivity of the AQD. 

In the high temperature regime ksT S> he/ a the phase coherence is reduced mainly 
due to thermal phonons. However, even at zero temperature coherence is lost due to 
purely quantum fiuctuations. This loss of coherence takes place on the same time scale 
as in the high temperature regime, but is incomplete in both 3D and 2D, which can be 
seen by evaluating the integral (j21|) in the limit t ^ a/c with T = 0. 



4. Measurement of the internal states 



The AQD dephasing can be detected in a Ramsey type experiment P^: The AQD is 
prepared in state |0) and hence initially decoupled from the BEC. A first 7r/2-pulse at 
t = changes the state to a superposition (|0) + |l))/-\/2 and a spin-echo-type vr-pulse 
at t = r/2 neutralizes the mean field shift. After a second 7r/2-pulse at t = r the AQD 
is found in state |1) with probability 

P(|l)) = i(l-e-^W). (26) 




Figure 3. The phase coherence e"''''^'^^ as a function of temperature T for interaction 
times Ti ~a/c (dotted), T2 = 2 cr/c (dashed) and T3 = lOcr/c (solid). In ID and 2D the 
interaction time r can be chosen to assure that dephasing changes significantly with 
temperature, whereas in 3D dephasing is independent of r for r 3> cr/c. The plot was 
produced with the same set of parameters as in figure ^ except for the temperature. 



However, this result is altered by decay of state |1) into state |0), atom loss, 
imperfect (noisy) detection of state |1), and additional dephasing due to environmental 
noise. We subsume these processes into three phenomenological constants, namely the 
detection probability P^, the probability of a spurious detection P^, and the dephasing 
rate 7^, which can all be determined experimentally. Taking these effects into account 
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we find an effective probability P to detect state |1) 

P(|l)) = iprf (1 - e-^^^)-^''^ + Ps (27) 
and the visibility V = {Pmax - Pmin)/{Pmax + Pmin) in the limit < 1 is given by 

where Pmax and Pmin correspond respectively to P(|l)) in the case of complete dephasing 
and no dephasing due to phase fluctuations. Kuhr et al recently demonstrated state- 
selective preparation and detection of the atomic hyperfine state for single cesium atoms 
stored in a red-detuned dipole trap. They showed that dephasing times of 146ms and 
ratios Ps/Pd of the order of 5 x 10~^ are achievable. If we choose 7^ = 10s~^, r = 10ms 
and Ps/Pd = 5 X 10~^ we find a visibility oiV = 69%, which shows that our scheme is 
feasible even in the presence of additional dephasing due to environmental noise. 



5. Conclusion 



We have shown that the dephasing of the internal states of an AQD coupled to a 
BEC is directly related to the temporal phase fluctuations of the BEC. Based on this 
relation we have suggested a scheme to probe BEC phase fluctuations nondestructively 
via measurements of the AQD coherences, for example, in a Ramsey type experiment. It 
was shown that the scheme works for a BEC with reasonable experimental parameters 
even in the presence of additional dephasing due to environmental noise. In particular, 
the scheme allows us to trace the dependence of the phase fluctuations on the trapping 
geometry of the BEC and to measure the BEC temperature. 

Our scheme is applicable even if the BEC is trapped in not strongly confined 
directions provided that the trapping potential is sufficiently shallow. We expect that the 
observed dephasing will be qualitatively different from our results only if the AQD size a 
is comparable to the classical harmonic oscillator amplitude ~ 1/u;. In addition, our 
findings indicate that the use of AQDs for quantum information processing, proposed 
recently in [221 12^] 5 may be constrained because of the unfavorable coherence properties 
of low dimensional BECs. 

The results in |17j suggest a natural extension of our work to entangled states 
between several AQDs, which might lead to a probe with higher sensitivity. However, the 
experimental requirements for the state-selective preparation and detection of entangled 
states are considerably higher than for a single AQD, which has to be considered in the 
analysis of an extended scheme. 
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